Stripe-tetragonal first order phase transition in ultra-thin magnetic films 
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We analyze the nature of the phase transition from a smectic stripe phase to a tetragonal phase 
predicted in analytic studies by Abanov et al. (Phys. Rev. B51, 1023 (1995)) and observed in 
experiments on ultra-thin magnetic films by Vaterlaus et al.(Phys. Rev. Lett. 84, 2247 (2000)). 
At variance with existent numerical evidence, for the first time we show results of Monte Carlo 
simulations on a two-dimensional model with competing exchange and dipolar interactions showing 
strong evidence that the transition is a weak first order one, in agreement with the theoretical 
predictions of Abanov et al. Besides the numerical evidence, we give further support to the first 
order nature of the transition analyzing a continuum version of the model and showing that it 
belongs to a large family of systems, or universality class, in which a first order transition driven by 
fluctuations is expected on quite general grounds. 
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With the advances in the experimental manipulation 
of materials at atomic length scales a renewed interest 
has grown in understanding the thermodynamic and me- 
chanical properties of systems such as ultra-thin films and 
quasi two-dimensional magnetic materials. Part of this 
interest is obviously motivated by the great amount of po- 
tential applications they find nowadays in many different 
technological fields (data storage, catalysis and electron- 
ics are only a few examples) . In this letter we are mainly 
concerned with the thermodynamic properties of ultra- 
thin magnetic films, like metal films on metal substrates 
(e.g. Fe on Cu @|, Co on Au @], see also Ref.[3] for a 
recent review on the topic). 

A large variety of magnetic films exhibit a spin re- 
orientation transition below some finite temperature Tr. 
That is, if the magnetic film is thin enough (a few atomic 
layers) the atomic magnetic moments tend to align in a 
direction perpendicular to the plane of the film, because 
the surface anisotropy overcomes the anisotropy of the 
dipolar interactions, which favors in-plane ordering. Un- 
der these circumstances the local magnetic moments can 
be regarded approximately as Ising variables. Any real- 
istic theoretical description of a magnetic thin film must 
include long-range dipolar interactions. The competition 
between exchange and dipolar interactions in these ma- 
terials give rise to stable modulated stripe-like patterns 
at low temperatures. In these states the magnetic mo- 
ments align along a particular axis forming ferromagnetic 
stripes of constant width h, so that moments in adjacent 
stripes are anti-aligned. Theoretical studies concluded 
that the stripe state in this systems is always the most 
stable one at low enough temperatures, provided that 
the exchangeparameter exceeds some small positive crit- 
ical valuejj, Calculations based on a continuum ap- 
proximation by Abanov et al. 0] predicted that, before 



reaching the paramagnetic state at high temperatures, 
the stripe phase undergoes a transition into a phase char- 
acterized by domains with predominantly square corners, 
that they called a tetragonal liquid. They also concluded 
that the stripe-tetragonal liquid transition should be ei- 
ther first order or the two phases might be separated by 
a third phase characterized by rotational domain walls 
defects, that they called an Ising nematic phase. Monte 
Carlo calculations on the square lattice carried out by 
Booth et al. 6] confirmed the presence of the tetragonal 
liquid phase, but they did not find any evidence of the 
Ising nematic phase. However, the results of Booth et al 
appeared to be consistent with a continuous transition 
rather that a first order one, as could be expected from 
Abanov et al theoretical results. Nevertheless, Booth et 
al pointed out that the possibility of a weak first order 
transition cannot be excluded on the base of their Monte 
Carlo simulations. Recent imaging studies using scan- 
ning electron microscopy with polarization analysis on 
ultra-thin films of fee Fe on Cu(100) verified the exis- 
tence of the tetragonal liquid phaseyl Q. No evidence 
was found of an intermediate nematic phase of the type 
predicted by Abanov et al.pj in the transition from the 
stripe phase to the tetragonal liquid. The thermody- 
namic nature of the stripe-tetragonal transition could not 
be determined by the imaging technique and this ques- 
tion remains unanswered. In this letter we present both 
Monte Carlo (MC) and analytical calculations that pro- 
vide new evidence that the stripe-tetragonal liquid transi- 
tion is indeed a weak fluctuation-induced first order one. 

We consider a system of magnetic dipoles on a square 
lattice in which the magnetic moments are oriented per- 
pendicular to the plane of the lattice, with both nearest- 
neighbor ferromagnetic exchange interactions and long- 
range dipole-dipole interactions between moments. The 
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FIG. 1: Energy per spin histograms for S = 2, L = 32 and 
different temperatures around the critical one T c w 0.79. The 
images below the histogram illustrate some typical equilib- 
rium spin configurations corresponding to the maxima and 
the minimum at the critical temperature. 



thermodynamics of this system is ruled by the dimen- 
sionless Hamiltonian [jj: 
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where S stands for the ratio between the exchange Jo > 
and the dipolar Jd > interactions parameters, i.e., 
S = Jo/Jd- The first sum runs over all pairs of near- 
est neighbor spins and the second one over all distinct 
pairs of spins of the lattice; ry is the distance, measured 
in crystal units, between sites i and j. The energy is mea- 
sured in units of Jd- The overall (known) features of the 
equilibriurnnhase diagram of this system can be found 
in Refs.p, IfiL Hc|. while several dyna mical propertie s 



at low temperatures can be found in R,efs. |llLll2lll.'ill< 
The threshold for the appearance of the stripe phase in 
this model is 5 C = 0.425 0,0. 

We performed extensive MC simulations of the model 
defined by Eq.flJ on a square lattice with N = L 2 
sites, for systems sizes ranging from L = 16 to 32 us- 
ing the Metropolis algorithm. Periodic boundary con- 
ditions were implemented using the Ewald summation 
technique. Although we obtained results for different 
values of S ranging from 1 to 3, most of the numerical 
work was focused on 5 = 2 (corresponding to a low tem- 
perature stripe phase of width h = 2), for which the 
first order nature of the transition is more clearly de- 
fined. To estimate both equilibration and decorrelation 
times we analyzed the behavior of the two-times corre- 
lation function for different system sizes and tempera- 
tures. After an equilibration period of up to 10 5 Monte 
Carlo Steps (MCS), every data point was calculated over 
a large single run for periods ranging from 2 x 10 6 MCS 
for the smallest size, up to 4 x 10 7 MCS for the largest 
one (L — 32). To locate the transition and characterize 
its order we calculated the energy per spin histograms for 
different temperatures, from which we obtained both the 
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order Binder cumulant[I3| V L = 1- (i? 4 ) /3 (H 2 )" . The 
typical behavior of the energy histograms is shown in 
FigQ]for 5 = 2 and L = 32. The double peak structure is 
characteristic of a first order phase transitional^]]. We 
observed that the measurement time (in MCS) needed to 
sample the energy distribution effectively increases ex- 
ponentially with the linear system size L. This is be- 
cause the system spends large periods of time in each 
one of the two phases, with a mean life time in ev- 
ery phase that increases with L, consistently with the 
expected[l^| increase in the energy free barrier between 
phases AF(L) ~ L . A careful inspection of the 
typical equilibrium spin configurations associated with 
every phase (i.e., those with energies around the max- 
ima of the energy distribution) shows indeed that the 
high temperature phase presents a tetragonal structure 
similar to that found by Booth et al. [6] (see Fig^l, 
while the low temperature phase is the h = 2 stripe one 
(also notice the coexistence of domains of both phases 
at energies corresponding to the minimum of the his- 
togram). Moreover, a MC calculation of the structure 



factor S(k) = (jJ2i a i e | ) at temperatures above, 
but near the critical one, show a very similar shape to 
that observed by Maclsaac et al. Q at higher values of 
5: four sharp peaks symmetrically placed on the two prin- 
cipal axes of the Brillouin zone, which characterize the 
tetragonal structure 0|. As the temperature increases 
the peaks become smeared into a circle with the shape 
of a four-peaked crown that gradually disappears. This 
indicates a continuous loss of the fourfold symmetry as 
the systems becomes paramagnetic. For large values of 
5 (6 > 3) this transition is also reflected in the presence 
of a secondary peak at high temperatures in the specific 
heat, which does not depend on the system size |fjj], in- 
dicating that the transition does not have an associated 
singularity in the thermodynamic potentials. As already 
observed by Booth et al |fj], this secondary peak becomes 
less pronounced as 5 decreases; we observed that for val- 
ues of S < 2.5 it becomes indistinguishable. However, 
the presence of the transition remains clearly detectable 
in the behavior of the structure factor. 

The temperature variation of the specific heat and 
the Binder cumulant for 5 = 2 and various system 
sizes is shown in Fig[2| We see that the location of 
the maximum of the specific heat and the minimum of 
the Binder cumulant shift in a size dependent fashion 
at pseudo-critical temperatures Tc X \l) and Tc 2 \l) re- 
spectively. Both quantities are plotted vs L~ 2 in FigG3 
showing the expected finite size scaling behavior for a 
temperature-driven first order phase transition^, 0] 
T C (1) (L) - T c + AL- d and T c (2) (£) - T c + BL~ d with 
B > A, where T c is the transition temperature of the 
infinite system. Note that the internal energies of both 
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phases (corresponding approximately to the energies of 
the maxima of the energy distribution) are located very 
close to each other. This property is also reflected in 
the rather shallow shape of the minimum of the Binder 
cumulant, evidencing the weak nature of the transition. 
We observed that these effects become more pronounced 
as S increases, with the internal energies of both phases 
approaching continuously to each other. For values of 
S > 2.6 the double peak structure of the energy his- 
togram (together with the minimum of the Binder cu- 
mulant) seems to disappear, or at least it becomes un- 
detectable for small system sizes, as can be appreciated 
in FigQ]for 5 — 2.6. We see that the internal energies 
of both phases (roughly corresponding to the energies 
of the maxima of the histogram) near the transition ap- 
proach to each other as 6 increases. This fact explains the 
seemingly continuous nature of the transition observed by 
Booth et al., whose calculations were performed for S > 3 
0. Clearly, considerably larger systems must be simu- 
lated in order to get reliable data for larger S. A similar 
effect is observed as S decreases approaching S c , where 
the double peak structure of the histogram also seems 
to disappear at all. Hence, the numerical data show an 
optimal value of S around 8 — 2, where the first order 
transition becomes strongest, that its, with the largest 
latent heat. Since strong finite size effects are always ex- 
pected in systems with long range interactions, numerical 
simulations in considerably larger systems are also re- 
quired in this case to get a more clear picture about this 
phenomenon. However, it is worth mentioning that some 
numerical simulations for larger system sizes (L = 48 and 



L = 64) showed that the double peak structure of the en- 
ergy histogram persists and becomes more pronounced as 
L increases. We now turn our attention to an analytic ap- 
proximation which gives us some insight in the expected 
outcome of such simulations for the other regions of the 
phase diagram. 

It was recognized long ago [lH . fill that systems in 
which the spectrum of fluctuations has a minimum at a 
non-zero wave vector can undergo a first order transition 
driven by fluctuations, in contrast to the second order 
transition predicted by mean field for this kind of sys- 
tems. Since the original work by Brazovskii, the proposed 
scenario was shown to describe correctly the phase transi- 
tions present in a large variety of systems like cholesteric 
liquid crystals, the nematic to smectic-C transition, pion 
condensates in neutron stars, onset of Rayleigh-Benard 
convection and micro-phase separation in diblock copoly- 
mers filj . HfSj . More recently the Brazovskii scenario 
has been successfully applied to the analysis of the phase 
transition between the disordered and modulated phases 
in three dimensional systems with attractive short range 
interactions and repulsive long range Coulomb interac- 
tions |2fl Elf . In spite of its success and ubiquity in cor- 
rectly describing the physics behind a phase transition in 
systems with competing interactions, the Brazovskii sce- 
nario was almost not considered in relation with ultra- 
thin films and dipolar systems • Indeed we have veri- 
fied that a continuum version of Hamiltonian presents 
a fluctuation induced first order phase transition for any 
value of 5. The starting point is a Landau-Ginzburg func- 
tional which in Fourier space has the form: 
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where the spectrum of fluctuations is given by A(k) — 
TQ+k 2 + J (k) / 5 and ro is proportional to the reduced tem- 
perature near the critical point of the mean field approx- 
imation. In our case the Fourier transform of the dipolar 
interaction is J{k) = X F 2 ({-1/2}, {1/2, 1}, -k 2 /4) - 
k, where i Fi ({a }, { /3i , /^j, z) is a generalized hyper- 
geometric series |23j. The spectrum of fluctuations 
presents a minimum at a nonzero k m i n — ko. This means 
that, for d > 2 the spectrum is minimized in a spherical 
shell in Fourier space. When fluctuations of the order pa- 
rameter are considered self-consistently, the degeneracy 
of the minimizing vector will make the disordered phase 
meta-stable for any finite temperature and drive the tran- 
sition to first order. More explicitly, within a Hartree ap- 
proximation the correlation function of the fluctuations 



in the disordered phase is given by: 

C-\k) = A(k) + 3u J j§^C(q) (3) 

Then, the renormalized mass r = r + 3u J C(q) is 
given by: 

f dq 1 
r = r Q + 3u / (4) 

J (2tt) 2 r + q 2 + J(q)/S v ' 

The point of absolute instability of the disordered solu- 
tion occurs at r = —k" 2 — J(ko)/5. It is easy to check that 
this point is only reached for ro — * — oo (i.e. for T — > 0), 
implying that the disordered phase never loses stability. 
At a finite temperature a modulated phase with a finite 
amplitude given by to 2 = (r m +kQ + J(ko)/S)/3u appears, 



4 



where the renormalized mass in the modulated phase r„ 
is given self-consistently by: 



= r + 3u 



+2 



k 



dq 



1 



(2tt) 2 r m + q* + J( ff )/5 
J(fco) 



(5) 



In a region r$ < Tq (S) a real solution to the above equa- 
tion exists. If a point Tq < rjj exists where the free ener- 
gies of the modulated and disordered solutions cross each 
other, then a first order transition driven by fluctuations 
appears. Following Brazovskii ^1 we have determined 
the point r§(5) where the free energy difference 
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FIG. 2: MC calculations for 5 = 2 and different system sizes 
I/, (a) Specific heat Cl vs. T; (b) Binder cumulant Vl vs. T. 



changes sign. It is possible to show that for any 5 > a 
first order transition appears. The transition line rg(<5) 
is a monotonously increasing function of S. 

In summary, we have shown strong numerical evidence 
that the phase transition between the low temperature 
stripe phase and the tetragonal phase observed in ultra- 
thin magnetic films is a weak first order transition. More- 
over, to the extent to which a continuum Hartree approx- 
imation of the model simulated is valid, the transition is 
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FIG. 3: Pseudo critical temperatures (maximum of the 
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specific heat) and T c (minimum of the Binder cumulant) vs. 



a first order one driven by fluctuations, at variance with 
the second order nature predicted by mean field theory. 
It would be interesting to perform simulations in larger 
lattices and with smarter Monte Carlo algorithms in or- 
der to analyze if the first order nature extends to a wider 
region of the phase diagram as implied by the continuum 
model. The challenge persists to probe the nature of this 
transition experimentally 
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FIG. 4: Energy per spin histograms for 8 — 2.6, L — 32 and 
different temperatures around the critical one. 
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